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[i]  Field  data  indicate  convincingly  that  the  drag  coefficient  of  the  ocean  surface  is  sea- 
state  dependent.  As  a  result,  under  steady  forcing  by  a  constant  wind  velocity  the  wind 
stress  on  the  ocean  surface  varies  with  time.  It  also  varies  with  space  if  the  wave 
development  is  limited  by  fetch.  A  quantitative  estimation  of  the  temporal  and  spatial 
variation  of  the  wind  stress  produced  by  a  constant  wind  velocity  is  presented.  The  method 
of  computation  combines  the  duration-  or  fetch-limited  growth  functions  of  wind¬ 
generated  waves  and  the  similarity  relation  of  the  ocean-surface  drag  coefficient  derived 
from  wavelength  scaling.  The  only  required  input  is  the  wind  speed.  The  results  indicate 
that  the  average  momentum  flux  from  atmosphere  to  ocean  is  much  larger  (about  50  to 
100  percent  higher,  and  especially  for  shorter  wind  events)  in  comparison  with 
calculations  using  the  drag  coefficient  or  dynamic  roughness  formulated  either  without  the 
wave  parameters  or  based  on  steady  state  wave  conditions. 
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1.  Introduction 

[2]  Over  the  last  several  decades,  many  wind  stress 
measurements  have  been  collected  from  oceans  and  lakes. 
These  measurements  are  used  to  investigate  the  dependence 
of  the  drag  coefficient  and  the  dynamic  roughness  of  the 
ocean  surface  on  various  environmental  parameters.  The 
reference  wind  speed  by  convention  is  U\o,  the  equivalent 
wind  speed  at  10-m  elevation  under  the  condition  of  neutral 
stratification.  In  the  intermediate  wind  speed  range  (approx¬ 
imately  from  7  to  20  m/s)  under  which  the  ocean  surface  is 
hydrodynamically  rough,  the  corresponding  drag  coeffi¬ 
cient,  C10,  displays  a  general  tendency  of  increasing  with 
wind  speed  [e.g.,  Garratt,  1977;  Wu,  1980],  Additional 
dependence  on  other  sea-state  parameters,  especially  the 
wave  age,  is  also  evident,  but  the  data  scatter  is  large  and  a 
quantitative  determination  of  the  additional  sea-state  depen¬ 
dence  remains  unsettled  [e.g.,  Donelan,  1990;  Toba  et  al., 
1990;  Geernaert,  1999;  Jones  and  Toba,  2001;  Drennan  et 
al.,  2003;  Smedman  et  al.,  2003;  Guan  and  Xie,  2004,  and 
references  therein].  Outside  the  intermediate  wind  speed 
range,  different  physics  govern  the  behavior  of  the  ocean- 
surface  drag.  For  lower  wind  speeds,  the  ocean  surface  is 
either  in  hydrodynamically  smooth  or  transitional  regime 
and  the  wave  influence  is  competing  with  the  viscous 
effects  [e.g.,  Donelan,  1990],  For  very  high  wind  speeds 
(say,  above  25  m/s),  extensive  breaking  occurs.  The  result¬ 
ing  spume,  flying  spray,  and  broad  regions  with  flow 
separation  act  like  a  shroud  shielding  the  fine-scale  wave 
roughness  from  the  airflow.  Field  measurements  indicate 
that  C10  reaches  a  maximum  near  35  m/s  and  then  decreases 
with  increasing  wind  speeds  [Powell  et  al.,  2003]. 
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[3]  The  following  discussions  focus  on  the  surface  wave 
effects  without  the  complication  of  extensive  wave  breaking 
that  leads  to  a  reduced  effective  surface  roughness,  and  limit 
the  scope  of  investigation  to  locally  generated  waves  by 
steady  winds  of  intermediate  velocities  and  neutral  stratifi¬ 
cation.  From  the  wave  dynamics  point  of  view,  the  more 
meaningful  reference  elevation  should  be  the  characteristic 
wavelength,  \p  [e.g.,  Kitaigorodskii,  1973;  Stewart,  1974; 
Donelan,  1990;  Makin  and  Kudryavtsev,  2002;  Oost  et  al., 
2002;  Hwang,  2004],  because  the  influence  of  surface 
waves  decays  exponentially  with  the  distance  from  the 
air-sea  interface  and  the  wavelength  serves  as  the  length 
scale  of  the  attenuation  rate  [e.g.,  Miles,  1957;  Phillips, 
1977],  Oost  et  al.  [2002]  presented  data  showing  an 
excellent  correlation  between  the  wind  friction  velocity, 
«*,  and  the  wind  speed  at  an  elevation  proportional  to  the 
wavelength,  UyJ2  or  U\,  for  measurements  under  neutral 
stratification  and  wind-sea  dominant  conditions;  the  corre¬ 
lation  coefficient  reaches  0.964.  Hwang  [2004]  carried  out 
an  analysis  of  the  drag  coefficient  Cv  2  referenced  to  Uyj2. 
The  wavelength  scaling  yields  considerable  improvement  in 
collapsing  field  measurements  of  the  drag  coefficient.  The 
assembled  data  set,  to  be  further  discussed  in  section  2, 
covers  a  wide  range  of  the  sea-state  conditions:  0.0235  < 
Wpujg  <  0.237,  0.0263  <  ujcp  <  0.237,  where  i0p  and  cp 
are,  respectively,  the  angular  frequency  and  the  phase 
velocity  of  the  wave  component  at  the  spectral  peak,  and 
g  is  the  gravitational  acceleration.  The  correlation  coeffi¬ 
cient,  Qc,  and  the  relative  root  mean  square  (RMS)  differ¬ 
ence,  S,  between  the  fitted  function  C\j2^puJg)  and 
measured  data  are  0.949  and  0.017,  respectively.  In  com¬ 
parison,  the  range  of  ( Qc ,  S)  for  various  C)0  functions  is 
(0.55-0.79,  0.025-0.034)  [Hwang,  2005a],  These  results 
strongly  suggest  that  the  similarity  relation  of  the  ocean- 
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Figure  1.  (a)  Drag  coefficient  C\j2(wpuJg)  measured  in  wind-sea  dominant  conditions  (DMAJ  data 
set),  and  dimensionless  roughness  (b)  k^0{^pujg)  and  (c)  z0/a(^pujg). 


surface  drag  coefficient  exists  in  wavelength  scaling.  For 
practical  applications,  it  is  necessary  to  convert  the  similar¬ 
ity  relation  of  (fy/2  to  C10,  which  is  not  a  trivial  problem. 

[4]  In  this  paper,  a  few  solutions  are  presented  to  quantify 
the  dependence  of  C10  on  wind  speed  and  sea  state  based  on 
the  similarity  relation  of  (fy/2  (section  2).  Because  the 
development  of  wind-generated  waves  is  duration-  and 
fetch-dependent,  it  follows  that  the  drag  coefficient  of  the 
ocean  surface  under  steady  forcing  by  a  constant  wind 
velocity  is  temporally  and  spatially  varying.  The  air-sea 
momentum  flux  computed  with  the  drag  coefficient  or 
dynamic  roughness  established  under  steady  state  wave 
conditions  or  without  considering  the  wave  factor  would 
underestimate  the  actual  magnitude  of  the  wind  input  into 
the  water  body.  Combining  the  similarity  relation  of  the 
drag  coefficient  derived  from  wavelength  scaling  and  the 
growth  functions  of  wind-generated  waves,  the  calculated 
C10  displays  significant  temporal  and  spatial  variations.  Its 
magnitude  is  about  a  factor  of  2  higher  in  the  first  couple  of 
hours  or  kilometers  compared  to  the  values  at  the  mature- 
wave  stage  (section  3).  Additional  discussions  on  the 
correlation  between  the  drag  coefficient  and  the  dynamic 
roughness  are  presented  in  section  4,  and  a  summary  is 
given  in  section  5. 

2.  Drag  Coefficient  of  the  Ocean  Surface 
2.1.  Wavelength  Scaling 

[5]  With  the  wavelength  scaling,  the  drag  coefficient  is 
represented  by  (fy/2  =  u\lU%j2.  Applying  this  scaling  to  field 
data  measured  under  the  condition  of  local  wind-wave 
generation  [Donelan,  1979;  Merzi  and  Graf,  1985;  Anctil 


and  Donelan,  1996;  Janssen,  1997]  (hereinafter  referred  to 
as  the  DMAJ  data  set;  the  experimental  conditions  were 
summarized  by  Hwang  [2004]),  the  data  scatter  of  (fy/2 
given  as  a  function  of  wpujg  is  markedly  reduced  in 
comparison  to  that  of  C]0  functions.  The  following  function, 
shown  as  the  solid  curve  in  Figure  la,  is  established  from 
the  DMAJ  data  set, 


with  Ac  =  1.220  x  10”2,  ac  =  0.704,  and  ( Qc ,  S)  =  (0.949, 
0.017). 

[6]  With  wavelength  scaling,  kpZ0  is  the  natural  expres¬ 
sion  of  the  dimensionless  roughness,  where  kp  is  the  wave 
number  at  the  spectral  peak.  Applying  the  function  of 
logarithmic  wind  speed  profile  to  equation  (1),  the  dimen¬ 
sionless  roughness  is  given  by 

kpz0  =  it  exp(-K,Cf°25J ,  (2) 

where  k,  =  0.4  is  the  von  Karman  constant.  Equation  (2)  is 
shown  as  the  solid  curve  in  Figure  lb.  For  comparison, 
the  dashed  curve  is  the  data-fitting  function  kpZ0  - 
1.889(uyVg)3'327,  with  (Qc,  S)  =  (0.949,  0.062). 

[7]  Figure  1  c  displays  the  result  of  normalizing  z0  with  the 
RMS  surface  displacement,  a.  The  dimensionless  roughness 
expressed  as  the  Chamock  parameter,  z0„  =  z0g/ul,  is 
more  scattered,  with  (Qc,  S)  =  (0.802, 0.159).  Because  z0/a  = 
kpZ0/s * 5,  where  s*  =  kpG2,  which  is  equivalent  to  e*(j* 
for  deep-water  waves,  e*  =  a2g2/(/j0,  and  u*  =  <s>pU\Jg, 
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Figure  2.  C i0  computed  from  equation  (1),  plotted  as  a  function  of  (a)  wpK*/g  and  (b)  wpU\cJg.  See  text 
for  the  description  of  the  Cv 2  curves. 


fetch-growth  functions  can  be  used  to  correlate  kpZ0  and  Zq/ct. 
Hwang  and  Wang  [2004,  2005]  developed  a  higher-order 
data-fitting  technique  to  obtain  the  empirical  power  law 
wave-growth  functions  with  variable  wave  development 
rates;  that  is,  the  exponents  of  the  power  law  functions 
may  vary  with  the  dimensionless  fetch,  duration,  or  wave 
frequency.  For  the  present  illustration,  the  growth  function  is 
represented  by  the  dependence  of  the  dimensionless  wave 
energy  on  the  dimensionless  wave  frequency, 

e*  =  Rw£.  (3) 

For  the  first-order-fitted  function,  R\  =  2.94  x  10~3,  rx  = 
-3.42,  and  for  the  second-order-fitted  function,  R2  = 
exp(-6.1384)w°'61021n“%  r2  =  -2.4019  -  1.22041nw* 
[Hwang,  2005b],  The  results  of  converting  kpZ0  to  zq/ct 
using  the  wave-growth  functions  are  shown  in  Figure  lc. 
For  comparison,  the  data-fitting  function  z0/cr  =  8.345(ujpm*/ 
g)2'863,  with  (Qc,  S)  =  (0.927,  0.099),  is  also  shown.  The 
results  illustrated  in  Figure  1  suggest  that  the  similarity 
relation  of  the  ocean-surface  drag  coefficient  exists  in 
wavelength  scaling,  and  that  equation  (1)  can  be  used  to 
represent  both  the  drag  coefficient  and  the  dynamic 
roughness  of  the  ocean  surface. 

2.2.  Conversion  to  10-m  Reference  Elevation 

[8]  For  practical  applications,  the  height  of  reference 
wind  speed  is  generally  set  at  10  m.  The  drag  coefficients 
C10  and  Cv 2  are  related  by 

C10  =  Cy2R2a,  (4a) 


R  ~Vm~ 


In  it  —  ln(fcpZo) 


In 


(  t4*gl0  \ 

\Cy/2kpz0U^  J 


(5) 


—  2  ln(Rj/) 


and  u>**  =  Upiijg.  It  is  of  interest  to  investigate  the  properties 
of  C10  as  a  function  of  wind  speed  and  sea-state  parameters 
based  on  the  similarity  function  of  C\j2  described  in  the  last 
section.  If  to**  and  kp  are  given,  the  procedure  to  derive  Cfo  is 
applying  equation  (1)  to  obtain  Cv2(uJ**),  then  equation  (2) 
to  obtain  ApZ0(w**),  and  finally  equation  (4b)  to  derive 
Cio(w**,  kp).  Of  more  interest  is  C10(u>**,  f710),  which  can  be 
derived  by  equations  (1),  (2),  (5),  and  (4a).  RLrcan  be  solved 
iteratively  for  given  to**  and  f/10.  From  numerical  experi¬ 
ment,  with  the  initial  guess  of  Ruo  =  1 ,  a  relative  error  of  1  %  is 
achieved  within  five  iterations.  Figure  2a  displays  the 
computed  CX0{upuJg)  for  a  range  of  w**  and  U\q.  For 
comparison,  the  C\j2(upuJg)  curve  (equation  (1))  is  also 
graphed  in  the  same  panel  with  a  thick  line  segment.  While 
C\/2  increases  monotonically  with  u)pw*/g,  C 10  curves  show 
local  maximum  near  u )pujg  =  0.25.  The  points  where  Ci0 
curves  intercept  the  C\j2  curve  correspond  to  the  condition 
that  the  peak  wavelength  is  20  m  and  C10  =  Cv 2-  Toward  the 
right-hand  side  (younger  waves)  of  the  interceptions,  C 10 
represents  the  drag  coefficient  referenced  to  a  wind  speed 
higher  than  Uyj2,  resulting  in  a  lower  magnitude.  Similarly, 
toward  the  left-hand  side  of  the  interceptions,  C/10  <  Uyj2 
therefore  C10  >  Cv2- 

[9]  Both  Cyj2(upuJg)  and  C'w(^pujg)  can  be  given  in 
terms  of  Wpt/io/g  by  the  substitution  of 


C 10 


NrT- 

k  \kpZa  J  \ 


(4b) 


WpU\D  _  WpU*/g 
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Figure  3.  (a)  Ci0  calculated  using  equation  (7)  and  presented  as  a  function  of  f/10;  (b)  same  as  Figure  3a, 
but  Ci0  is  calculated  using  equation  (1).  The  DMAJ  data  are  superimposed  for  comparison. 


The  resulting  C10(Wpt/10/g)  curves  for  Ul0  =  5,  10,  15,  and 
20  m/s  are  shown  in  Figure  2b.  Also  superimposed  in  this 
figure  are  the  corresponding  Cyj2(upU]0/g)  curves,  shown  as 
dotted  lines  that  increase  monotonically  with  u)pUu)/g.  The 
thick  line  segment  plotted  in  Figure  2b  is  the  empirical 
function  derived  from  data  fitting  of  the  DMAJ  data  set 

(7, 

with  A io=  1.289  x  10”3,  ai0  =  0.815,  and  (Qc,  S)  =  (0.893, 
0.024)  [Hwang,  2005a], 

[10]  From  the  point-of-view  of  practical  applications, 
equation  (7)  is  more  convenient  than  equation  (1)  for 
obtaining  Ci0  using  the  similarity  relation  of  the  drag 
coefficient  derived  from  wavelength  scaling.  If  gj*  and  kp 
are  given,  the  procedure  to  obtain  Ci0(go*,  kp)  is  straightfor¬ 
ward,  through  equations  (7),  then  (2)  and  (4b).  To  derive 
Cio(w*,  U\0),  only  a  slight  modification  of  the  procedure  is 
needed:  kp  can  be  calculated  from  the  dispersion  relation 
because  gj/;(=  oj*g/f/i0)  is  available.  Also,  C]0( to*,  U\ 0)  can 
be  constructed  from  equation  (7)  with  wp  and  U\0  input.  For 
intermediate-  and  shallow-water  wave  conditions,  the  addi¬ 
tional  information  of  water  depth  is  needed  to  obtain  kp  from 
wp  to  apply  equation  (4b).  Numerical  computations  show 
that  equations  (1)  and  (7)  produce  very  similar  C10  at  about 
10  m/s  wind  speed.  At  higher  wind  speeds,  equation  (7) 
yields  lower  C10  compared  to  equation  (1),  and  vise  versa  at 
lower  wind  speeds.  The  difference  is  generally  within  25% 
of  each  other.  Figure  3a  shows  C\0  curves  computed  with 
equations  (7),  (2),  and  (4b).  Superimposed  on  these  curves 
are  the  DMAJ  data  sorted  into  three  gj*  bins  and  ensemble- 
averaged  over  10  wind  speed  sub-ranges  for  each  gj*  bin  to 
match  the  computational  curves.  For  comparison,  the  results 
computed  using  equations  (1),  (2),  (5),  and  (4a)  are  displayed 
in  Figure  3b  with  DMAJ  data  superimposed. 


[n]  In  the  literature,  C,0  is  frequently  expressed  as  a 
linear  function  of  £/10.  The  slope  of  the  linear  function  is 
generally  steeper  for  younger  seas  and  milder  for  more 
mature  seas  [e.g.,  Guan  and  Xie,  2004],  a  feature  that  is 
consistent  with  the  result  shown  in  Figure  3.  It  should  be 
pointed  out  that  representing  Ci0  as  a  linear  function  of  Ui0 
is  dimensionally  inconsistent,  and  it  yields  a  constant  wind 
stress  for  a  constant  wind  speed.  Although  the  sea-state 
influence  can  be  included  through  introducing  sea-state 
dependence  for  the  slope  of  the  linear  function,  it  remains 
unsatisfactory  that  such  an  expression  is  dimensionally 
inconsistent. 


3.  Temporal  and  Spatial  Variation  of  the  Drag 
Coefficient 

[12]  Because  wave  development  is  dependent  on  fetch 
and  duration,  under  steady  wind  forcing  the  wind  stress 
exerted  on  the  water  surface  varies  with  time  and  space  as  a 
consequence  of  the  sea-state  dependence  of  the  drag  coef¬ 
ficient.  The  spatial  and  temporal  variation  of  the  wind  stress 
can  be  quantified  through  the  fetch-  and  duration-limited 
wave-growth  functions 


(  Axi 

\  Ptg 


where  x*  =  xg!U\o,  x  is  fetch,  t*  =  tg!Uw,  and  t  is  duration. 
The  coefficients  of  the  duration-limited  growth  functions 
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Figure  4.  Spatial  and  temporal  variation  of  Cx0  under  steady  wind  forcing  in  (a,  c)  dimensionless  fetch 
and  duration,  respectively,  and  (b,  d)  in  dimensional  fetch  and  duration,  respectively. 


can  be  written  in  terms  of  the  coefficients  of  the  fetch-  can  be  calculated  with  equations  (7),  (2),  and  (4b).  Figure  4 
limited  growth  functions  [Hwang  and  Wang,  2004,  2005],  shows  Ci0  as  a  function  of  dimensionless  or  dimensional 

fetch  and  duration  for  wind  speeds  between  5  and  20  m/s.  A 
[*,(&  + 1)1*  a  much  larger  drag  coefficient  at  the  early  stage  of  wave 

=  A  __  p  =  _  development  is  evident.  The  average  wind  stress  in  the  first 

_i_  b  (9)  2  hours  of  a  wind  event  is  about  50  to  70%  higher  than  that 

Q  =  [ssRc(6  +  1)1*"  q  =  — — - ,  at  the  mature-sea  stage. 


where  Rc  is  the  ratio  between  the  effective  group  and  phase 
velocities  of  the  wave  component  at  the  spectral  peak.  For 
wind  seas,  Rc  ss  0.4  [Yefimov  and  Babanin,  1991], 
Generally,  the  fetch-  and  duration-limited  wave  growth 
data  can  be  represented  by  first-order-fitted  power  law 
functions  (with  constant  exponents)  [Hwang  and  Wang, 
2004,  2005]  for  the  fetch  range  ~102  <  jc*  <  ~104. 
The  coefficients  are  Ax  =  6.19  x  10-7,  ax  =  0.811,  Bx  = 
11.86,  bx  =  -0.237,  Px  =  1.27  x  10~8,  /?,  =  1.062,  Qx  = 
36.92,  and  qx  =  —0.310.  For  a  broader  coverage  of  x*  the 
development  rate  is  not  constant.  The  variable  development 
rate  can  be  obtained  from  data  by  using  a  higher-order 
polynomial  function  in  log-log  scales  for  data  fitting.  To  the 
second  order,  the  coefficients  of  the  fetch  laws  become 

A2  =  e°«x*  ~  a2  =  OL\  +  2ol2\wc* 

a  ,  (10) 

B2=e?°xT-  *  =  Pi  +  2(32lnx„, 

with  c«o  =  3.0377,  ai  =  —0.3990,  a 2  =  0.0110,  (30  = 
—  17.6158,  Pi  =  1.7645,  and  p2  =  —0.0647.  The  correspond¬ 
ing  coefficients  of  the  second-order-fitted  power  law 
functions  for  e*(t*)  and  w*(f*)  can  be  calculated  by 
equations  (10)  and  (9).  Applying  the  above  wave-growth 
functions  to  derive  the  dimensionless  frequency,  the  spatial 
and  temporal  evolution  of  C10  under  steady  wind  forcing 


4.  Discussion 

[13]  There  has  been  considerable  uncertainty  about  the 
quantitative  properties  of  the  drag  coefficient  of  the  ocean 
surface.  It  was  initially  considered  to  be  a  constant  and  later 
revised  to  be  increasing  with  wind  speed.  Although  it 
remains  unsettled  on  how  the  drag  coefficient  increases 
with  wind  speed,  field  data  suggest  that  the  rate  of  increase 
is  sea-state  dependent.  Another  approach  to  studying  the 
sea-surface  drag  is  through  the  investigation  of  the 
dynamic  roughness.  (This  seems  to  be  a  unique  feature  of 
the  ocean-surface  drag;  almost  all  other  branches  of  fluid 
mechanics  describe  the  surface  drag  by  the  drag  coefficient, 
expressed  correctly  in  dimensionally  consistent  functions.) 
The  Chamock  parameter,  z0*  =  z0g/w*,  was  originally 
considered  to  be  a  constant  [Charnock,  1955],  but  later 
was  also  found  to  be  sea-state  dependent.  Donelan  [1990] 
normalized  z0  by  a.  The  resulting  function  is  much  more 
successful  in  explaining  measurements  from  disparate 
experiments.  While  there  is  good  progress  in  characterizing 
z0,  it  remains  a  puzzle  that  the  behavior  of  its  counterpart, 
the  drag  coefficient,  is  so  difficult  to  understand.  The 
dynamic  roughness  z0  is  not  a  measured  quantity.  Its 
derivation  depends  on  the  application  of  the  logarithmic 
wind  profile,  which  can  be  written  as 

cr0  5=-  In  —  ,  (11) 
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where  Cz  is  the  drag  coefficient  referenced  to  the  wind 
speed  at  elevation  z.  The  correlation  of  the  drag  coefficient 
and  the  dynamic  roughness  should  be  one-to-one  and 
deterministic  as  elicited  by  equation  (11). 

[14]  Hwang  [2004]  suggested  that  the  difficulty  in  finding 
the  similarity  properties  of  the  ocean-surface  drag  coeffi¬ 
cient  can  be  attributed  to  the  choice  of  the  arbitrary  10  m  as 
the  length  scale  for  wind  speed  reference.  Indeed,  when 
processed  with  wavelength  scaling,  the  wind  stress  mea¬ 
surements  in  oceans  or  lakes  under  wind-sea  dominant 
conditions,  either  represented  by  the  drag  coefficient  or  by 
the  dynamic  roughness,  show  strong  similarity  behavior  on 
upujg  (Figure  1). 

5.  Summary 

[15]  The  drag  coefficient  is  sea-state  dependent,  so  al¬ 
though  the  forcing  wind  velocity  is  constant  and  steady,  the 
wind  stress  on  the  ocean  surface  varies  with  time  and  space. 
The  temporal  and  spatial  variation  of  the  wind  stress  repre¬ 
sented  by  C10  is  quantified  by  combining  the  similarity 
relation  of  the  drag  coefficient  and  the  growth  functions  of 
wind-generated  waves.  The  results  show  that  the  wind  stress 
is  much  larger  in  the  first  few  hours  of  a  wind  event  than  that 
at  a  later  time  when  the  wave  field  is  more  mature.  Wind 
stress  computations  using  the  drag  coefficient  or  dynamic 
roughness  formulated  without  considering  the  wave  factors 
or  based  on  steady  state  wave  conditions  may  underestimate 
the  wind  input  into  the  ocean  by  as  much  as  a  factor  of  2 
(Figure  4).  The  procedure  described  in  this  paper  provides  a 
parametric  means  to  estimating  the  wave  influence  on  the 
momentum  transfer  in  a  coupled  air-sea  interaction  system  in 
the  absence  of  a  numerical  ocean  wave  model.  It  is  relatively 
easy  to  use  and  requires  only  the  wind  speed  input.  The  sea- 
state  information  is  embedded  in  time  and  space  through  the 
wave-growth  functions.  The  procedure  outlined  in  this  paper 
is  applicable  to  wind-generated  wave  conditions  in  the  field 
environment.  The  drag  coefficient  and  dynamic  roughness 
under  mixed-sea  conditions  remain  difficult  to  parameterize 
at  this  stage. 
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